METHOD FOR REMOVAL OF PID DYNAMICS FROM MPC MODELS 

************************* 
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Not applicable 
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BACKGROUND OF THE INVENTION 

[0001] Model Predictive Control (MPC) refers to a class of algorithms that compute a sequence of 
manipulated variable adjustments in order to optimize the future behavior of complex multivariable 
processes. Originally developed to meet the needs of petroleum refineries and chemical processes, 
MPC can now be found in a wide variety of application areas including chemicals, food processing, 
automotive, aerospace, metallurgy, and pulp and paper. A well-known implementation of MPC in 
chemical and refinery applications is Dynamic Matrix Control or DMC. 

[0002] The MPC Controller employs a software model of the process to predict the effect of past 
changes of manipulated variable and measurable disturbances on the output variables of interest. 
The independent variables are computed so as to optimize future system behavior over a time 
interval known as the prediction horizon. In the general case any desired objective function can be 
used for the optimization. The system dynamics are described by an explicit process model, which 
can take, m principle, a number of different mathematical forms. Process input and output 
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constraints are included directly in the problem formulation so that future constraint violations are 
anticipated and prevented. 

[0003] In practice a number of different approaches have been developed and commercialized in 
implementing MPC Controllers. The most successful implementations have made use of a linear 
model for the plant dynamics. The linear model is developed in a first step by gathering data on the 
process by introducing test disturbances on the independent variables and measuring the effects of 
the disturbances on the dependent variables. This initial step is referred to as Identification and the 
novel use of this identification data is the essence of this invention. 

[0004] U.S. Patents 4,349,869 and 4,616,308 describe an implementation of MPC control called 
Dynamic Matrix Control (DMC). These patents describe the MPC algorithms based on linear 
models of a plant and describe how process constraints are included in the problem formulation. 
Initial identification of the MPC controller using process data is also described. 

[0005] By way of further background this Identification of process dynamics requires a pre-test in 
which the independent variables of the process are moved in some pattern to determine the effect 
on the dependent variables. In a chemical or refinery process the independent variables include 
the PID (proportionaWntegral-derivative) controller set points for selected dependent variables, the 
valve positions of PID controllers in manual, and temperatures, material flows, pressures and 
compositions that are determined outside the scope of the controller's domain. For any process 
Identification test, the independent variables are fixed for the analysis of the data. Further the tuning 
of any of the PID controllers in the domain of the MPC controller is fixed. The MPC controller that is 
built to use the dynamic process models from the Identification must have exactly the same 
configuration of independent variables that existed when the Identification was performed. Thus the 
PID controller configuration that is present during Identification imbeds the PID controller dynamics 
in the dynamic model. 

[0006] This characteristic of current Identification technology represents an unsolved problem that 
is addressed by this invention. The problem creates a limitation on the use of MPC technology that 
manifests itself in two different areas. 

[0007] The first application area is MPC itself. Because the dynamics of the PID controllers are 
imbedded in the MPC model, any change in the tuning of a PID controller or changing of the PID 
state from auto to manual or vice versa changes the dynamic model. To correct this it has been 
required to retest the process unit with the changed conditions. A well-designed Identification test 
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for a complex multivariable process will be a 2-3 week effort with skilled people, costing two to three 
hundred thousand dollars. 

[0008] The second application area is in the field of Operator Training Simulators. Effective training 
simulators are important to the chemical process industry. The large investments in new chemical 
processes and the safety implications of the complex processes require a well-trained operator 
group. This is important especially for process units that remain on computer control for extended 
periods of time, since the operators do not have the opportunity to control the unit. MPC models 
have not been used in creating training simulators because of the aforementioned problem that the 
PID controller configuration that is present during Identification imbeds the RID controller dynamics 
in the dynamic model. The result of this is that authentic training is difficult because the operators 
cannot change the state (auto or manual) of the PID controllers without reducing the fidelity of the 
model. Current state of the art training simulators are based on approximated first principle 
engineering equations that have been simplified to such a degree that they do not represent the 
actual process dynamics. A survey of control rooms in the chemical process industry will reveal that 
they are rarely used after a stait-up as the operating personnel learn that the simulator does not 
reflect the actual dynamic behavior. A training simulator based on an identification model that has 
the fidelity to hold a process at constraints, display all temperatures, pressures, flows, and valve 
positions and allow the operator to switch any PID controller to manual or auto would be a powerful 
tool for training. 

[0009] Numerous unsuccessful attempts have been made by practitioners in the field to address 
this identification issue. One approach would be to run the identification test with the regulatory 
control scheme in manual. This of course fails because there are not enough operators available to 
handle a complex process with 50 to 100 PID controllers on manual. Other attempts have been 
made to conduct a standard identification test with the regulatory scheme in place but to then set up 
the model with the valve positions as the independent variables. These approaches always lead to 
failure with erratic results. It has become recognized that this approach fails because the valve 
positions are correlated through the PID dynamics via measured and unmeasured disturbances that 
are always present in a real world identification test and are thus not independent. 

[00010] The recognition of this fact and the method of removing the noise and unmeasured 
disturbances from the data set are the essence of this invention. The method of removing the noise 
and unmeasured disturbances is a procedural one. The process is first identified with a commercial 
tool such as AspenTech's DMC-Model. This step breaks the correlation between the unmeasured 
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disturbances, the noise, and the PID controllers. The PID controller dynamics can then be removed 
from the resulting model without interference from the noise and unmeasured disturbances. 



BRIEF SUMMARY OF THE INVENTION 

[00011] An object of this invention is to provide a method for removing the dynamics of the PID 
controllers from the MPC controller used in a multivariable control process. This allows creation of 
valve-based off-line process simulators. 

[0001 23 It is a further object of this invention to provide such a method that can be used in various 
implementations of MPC controllers. 

[00013] It is a further object of this invention to provide a method to generate new MPC controllers 
for complex multivariable process control when a change has been made in any PID control 
configuration or tuning and to do so without having to conduct new identification testing of the 
process. 

[00014] It is a further object of this invention to generate a process simulator based on valve 
positions with the effect of unmeasured disturbances removed so that a high fidelity process 
simulator Is available for process simulation and training. Such a simulator could be used for 
simulation in any controller configuration and with various tuning configurations on each individual 
controller. 

[00015] In accordance with this invention there is provided a method used in model predictive control 
applications for removing PID controller dynamics from a controller model of a process having a 
plurality of independently controllable, manipulated variables and at least one controlled variable 
dependent upon the independently controllable, manipulated variables that includes at least the 
steps of: gathering data about the process by separately introducing a test disturbance in each of 
the manipulated variables and measuring the effect of the disturbances on the controlled variable; 
using the effects of the disturbances on the controlled variable to generate a first linearized matrix 
model relating the at least one controlled variable to the independently controllable, manipulated 
variables; interchanging selected valve position controlled variables with their corresponding 
selected independently controllable, manipulated PID controller set point variables in the linearized 
model using matrix row elimination mathematics to generate a second linearized model that has a 
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new set of independently controllable, manipulated variables, the second model having the 
dynamics of the selected independently controllable, manipulated PID controller set point variables 
removed from the model. 

[00016] To use this model in a control context the method includes controlling a process having a 
plurality of independently controllable, manipulated variables and at least one controlled variable 
dependent upon the independently controllable, manipulated variables including the steps of: 
gathering data about the process by separately introducing a test disturbance in each of the 
manipulated variables and measuring the effect of the disturbances on the controlled variable; using 
the effects of the disturbances on the controlled variable to generate a first linearized dynamic 
model relating the at least one controlled variable to the independently controllable, manipulated 
variables; interchanging selected valve position controlled variables with their corresponding 
selected independently controllable, manipulated PID controller set point variables in the first 
linearized dynamic model using matrix row elimination mathematics to generate a second linearized 
dynamic model that has a new set of independently controllable, manipulated variables, the second 
linearized dynamic model having the dynamics of the selected independently controllable, 
manipulated PID controller set point variables removed from the second linearized dynamic model; 
measuring the present value of the variables; calculating for discrete intervals of time from the 
gathered data about the process, the measured present values and pre-selected operating 
constraints a set of moves for present and future times for at least the manipulated variables to 
obtain new values for the manipulated variables and to move the at least one dependent 
controllable variable towards at least one of the constraints; and changing the process by adjusting 
the manipulated variables for the set of moves for present and future times to cause the process to 
move the at least one dependent controllable variable towards at least one of the constraints. 

[00017] To use this invention to generate new MPC controllers for complex multivariable process 
control when a change has been made in any PID control configuration or tuning and to do so 
without having to conduct new identification testing of the process the following process can be 
used: by interchanging at least one PID controller set point variable in an original linearized 
dynamic model with its corresponding valve position controlled variable in the original linearized 
dynamic model using matrix row elimination mathematics to generate a secondary linearized 
dynamic model that has the at least one corresponding valve position as a new independently 
controllable, manipulated variable; then externally emulating new desired PID tuning via 
mathematical emulator to emulate the effect of the at least one PID controllers new tuning with the 
secondary linerarized dynamic model; then testing the secondary linearized dynamic model with it's 
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emulated PID tuning by stepping each of it's manipulated variables to obtain the new linearized 
dynamic model that will now contain the dynamics of the at least one PID controllers. 

[00018] It should be noted that a regulatory control scheme can be easily emulated external to the 
process model via a DCS console or console emulator available in modern control packages. This 
allows the operator to put PID controllers in Manual-mode, break cascades, retune PID controller, 
or even re-configure the regulatory control scheme. 

[00019] To use this invention to generate a process simulator based on valve positions with the 
effect of unmeasured disturbances removed so that a high fidelity process simulator is available for 
process simulation and training the following method is used: first gathering data about the process 
by separately introducing a test disturbance in each of the manipulated variables and measuring the 
effect of the disturbances on the controlled variable; then using the effects of the disturbances on 
the controlled variable to generate a first linearized dynamic model relating the at least one 
controlled variable to the independently controllable, manipulated variables; then interchanging 
each independently controllable, manipulated PID controller set point variable with its 
corresponding valve position controlled variable in the first linearized dynamic model using matrix 
row elimination mathematics to generate a second linearized dynamic model that has the 
corresponding valve positions as a new set of independently controllable, manipulated variables, 
the second linearized dynamic model having the dynamics of the selected independently 
controllable, manipulated PID controller set point variables removed from the second linearized 
dynamic model; then externally emulating a desired regulatory control schemes via mathematical 
emulators to emulate PID controllers in either manual, cascade, or automatic modes. As before it 
should be noted that a regulatory control scheme can be easily emulated external to the process 
model via a DCS console or console emulator available in modern control packages. This allows 
the operator to put PID controllers in Manual-mode, break cascades, retune PID controller, or even 
re-configure the regulatory control scheme 

[000201 The most common method of Identification currently used in oil refining and chemical 
processes is the Dynamic Matrix Identification (DM I). DM I will be used to illustrate the methodology 
of this invention, but it should be understood that the invention is not limited to a specific 
Identification technique. 
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DESCRIPTION OF DRAWINGS 



[00021] Figure 1 is a schematic demonstrating dependent and independent variables 

[00022] Figure 2 is a more complex flow schematic of a fractionator 

[00023] Figure 3 is a simulation of the fractionator model based on valve positions 

[00024] Figure 4 demonstrates the results from a plant test of the fractionator 

[00025] Figure 5 is a simulation of the fractionator with the PID controllers 

[00026] Figure 6 is a demonstration of the fractionator with the original and recovered values 



' DETAILED DESCRIPTION OF INVENTION 



[00027] The invention is a method used in conjunction with model predictive control for removing the 
dynamics of PID controllers from MPC controllers. 

[00028] An MPC process model is a set of linear equations so it should be mathematically possible 
to interchange any independent variable with a dependent variable provided a relation exists 
between the independent and dependent variable. A candidate set for that transformation is the set 
point (independent) for a PID controller and the associated valve position (dependent) for that PID 
controller. 

[000291 An MPC controller is often based on a linear model of a process system. For example 
see Figure 1 . Although the invention to be described here has applications in many fields the 
examples used will be from chemical and refinery process applications. 

[00030] There are two types of variables in any system; the independent variables and the 
dependent variables. The independent variables are inputs to the system. The independent 
variables are further divided into manipulated and disturbance (feedforward) variables. Manipulated 
variables are those that can be changed by the human operator, such as valve positions or PID 
controller set points. The furnace fuel gas flow set point (1) is a manipulated independent variable. 
Disturbance variables are those independent variables that have an effect on the system, but 
cannot be changed by the human operator. Variables such as feed composition, feed temperature, 
and ambient temperature are examples of disturbance variables. The furnace inlet temperature (2) 
is a disturbance variable. 
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[00031] Dependent variables are outputs from the system. Dependent variables are affected by 
changes in the independent variables. The human operator cannot directly change them. The 
values of dependent variables can be controlled, however, by correctly changing the values of the 
manipulated variables. Further, as disturbances enter the system, the manipulated variables must 
be correctly adjusted to counteract the disturbance. The furnace outlet temperature (3) and stack 
temperature (4) are examples of dependent variables. 

[00032] The use of linear models allows the use of matrix mathematics in describing complex and 
multivariable control. There are several general formulations of MPC models. A general model for 
control is the step response model: 

SO, =A,^,M, +-'-hA,jAIj +'-+4.««^a7„„,, 

SO. -4iA/i +.-+4-^-A/^ +-+A.nin^^n.d 

Equation 1 : Step Response Dynamic Matrix, Block Matrix Form 



where, 



SO.= 





-0,0 




-0,0 




-0,0 



, the accumulative change in the dependent variable at each time step, 



AIj = 



AI 



AI 



AI 



y.3 



AI. r 



, the step change in the/'' independent variable at each time step, and 



a 



J^tjMCoef ^iJXncoef-X) ^ij\incoef-2) 



, the Dynamic Matrix. 
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[00033] An alternate form of this Step Response equation is the Finite Impulse Response (FIR) form. 
It can be derived from the Step Response form as described below. 

[00034] Recalling from the definitions that: 

hjjc=<^ij,k ior k = \. 

Km = - ^i,Mk-i) foTk:2^ ncoef 
and that 

^O^jc = - O,. for A; : 1 ncoef 

we can difference the above system of equations to give: 

AO, = B,,AI, - +B,jAl - +B,^^Ai„^ 

AO, = B,,AI, - +B^AIj - +B,„^Ai„^^ 
^O^p B^paA "■+ Br^j^ij -+ B^^„,^AI^ 



Equation 2: Finite impulse Response Equations - Block Matrix Form 
where 





o,,-o,^. 










A0,= 


0■,^-0,^ 


, the change in the dependent variable across each time interval, 




_Pi,ncoef ~ ^ilncoefA) _ 





AIj 



AI 



A/ 

AI 



ja 



y.3 



j, ncoef 



as above, and 
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Kj,2 



Km 



, the model matrix of Impulse Coefficients. 



[00035] There are five forms of these equations, and we have shown only the first two. While these 
forms are mathematically equivalent, and while all forms may be used for identification prediction 
and control, they have very different properties. 

SO- AM - Most often used for control calculations. 
O AO = BM - Used for identification of steady state variables. 
M AAO = BAAI - Used for identification of ramp variables. 
3 so = BSI - Not commonly used. Old IDCOM control formulation. 

AO = AAAI - Not commonly used. 

^ [00036] C. R. Cutler and C. R. Johnston discuss the properties of these forms of the matrix in a 
I'* paper, "Analysis of the Forms of the Dynamic IVIatrix", in the Proceedings of the Instmment Society 
% of America ISA 85 Advances in Instrumentation Volume 40, Number 1 - October 1985. 

[00037] The use of these linear modeling techniques, including the identification of the model and 
the use of the model for control and the use in control with constraints is described in two U.S 
patents, 4,349,869 and 4,616,308. These patents are incorporated herein by reference. 

[00038] We will now derive the algorithm of this invention to demonstrate the removal of the PID 
dynamics from the controller. The derivation is from the FIR model of equation 2. To derive the 

algorithm, we assume that the j independent variable is the set point of a PID controller and the 

i^^ dependent is the PID valve response to that set point change. We wish to re-constitute the 
model so that the valve is the independent variable in the process model; that is to say, we wish to 
remove the dynamics of this PID controller from all affected model responses. This can be 

accomplished by interchanging the dependent variable with the j^^ dependent variable, as 

follows: 
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the Identity Matrix. 



Note that this is nothing more than equation 2 above with an Identity matrix multiplying theAO' j 
By performing row elimination operations (pivoting), we get; 
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Which can be re-written as: 
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Which can be rearranged to; 
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or reassembling the matrix equation we get; 
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[00039] Note that AO. and A/, have been interchanged so that the valve position is now an 

independent variable and the PID set point is now a dependent variable. This illustrates removing 
the PID dynamics from only one PID controller, but the algorithm is clearly general in that multiple 
independent/dependent variable pairs can be interchanged to remove the dynamics for multiple 
controllers. 

[00040] By way of further illustration a numerical example problem will now be illustrated to show 
how this approach is applied to a model predictive controller to remove the dynamics of a particular 
PID controller. 

[00041] Given an FIR model with two (2) independent variables, two (2) dependent variables and 
four (4) model coefficients, where the second independent variable is the set point of a PID 
controller and the second dependent variable is the valve position of the PID controller, we wish to 
re-constitute the model with the PID valve position as an independent variable instead of the PID 
set point. This requires that the dynamics of the PID controller be removed from all system 
responses according to the algorithm previously discussed. This example is also valid for the 
AO == BM^ , SO = BS I ,and AAO = BAAI forms of the equation. 



Dependent Var-1 

Independent Var-1 Independent Var-2 



^1,2,1 " 


0.5 


^,2,2 ~ 


0.4 


^1,2,3 ~ 


0.2 


^1.2,4 ~ 


0.1 
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Dependent Var-2 
Independent Var-1 

hxr -0-3 



Independent Var-2 

K2,= 0-25 

^2.2.4= 0-05 



[00042] The problem is specified in tiie matrix below. 
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Multiply Equation-5 by (-1/0.75) 
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Multiply Equation-5 by 0.5, add it to Equation-1 and replace Equation-1 
Multiply Equation-5 by 0.4, add it to Equation-2 and replace Equation-2 
Multiply Equation-5 by 0.2, add it to Equation-3 and replace Equation-3 
Multiply Equation-5 by 0.1 , add it to Equation-4 and replace Equation-4 
Multiply Equation-5 by 0.25, add it to Equation-6 and replace Equation-6 
Multiply Equation-5 by 0.15, add it to Equation-7 and replace Equation-7 
Multiply Equation-5 by 0.05. add it to Equation-8 and replace Equation-8 
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0 


0.15 


0.25 


0.75 


0 


0 


0 


0 


-0.067 


0 


0 


1 


Multiply Equation-6 by (-1/0.75) 


1.7 


0 


0 0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 


0 


0 


0 


0.76 


1.5 


0 0 


0 


0.5 


0 


0 


0 


1 


0 


0 


-0.533 


0 


0 


0 


0.28 


0.6 


1.5 0 


0 


0.4 


0.5 


0 


0 


0 


1 


0 


-0.267 


0 


0 


0 


0.14 


0.2 


0.6 1.5 


0 


0.2 


0.4 


0.5 


0 


0 


0 


1 


-0.133 


0 


0 


0 


0.4 


0 


0 0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1 .333 


0 


0 


0 


0.4 


0.4 


0 0 


0 




0 


0 


0 


0 


0 


0 


0.444 


-1 .333 


0 


0 


-0.04 


-0.4 


-0.3 0 


0 


0.25 


0,75 


0 


0 


0 


0 


0 


-0.2 


0 


1 


0 


-0.03 


-0.1 


-0.4 -0.3 


0 


0.15 


0.25 


0.75 


0 


0 


0 


0 


-0.067 


0 


0 


1 



Multiply Equation-5 by 0.5, add it to Equation-2 and replace Equation-2 

Multiply Equation-5 by 0.4, add it to Equation-3 and replace Equation-3 

Multiply Equation-5 by 0.2, add it to Equation-4 and replace Equation-4 

Multiply Equation-5 by 0.25, add it to Equation-7 and replace Equation-7 

Multiply Equation-5 by 0.15, add it to Equation-8 and replace Equation-8 



1.7 


0 


0 0 


0 


0 


0 




1 


0 


0 


0 


-0.667 0 


0 


0 


0.96 


1.7 


0 0 


0 


0 


0 




0 


1 


0 


0 


-0.311 -0.667 


0 


0 


0.44 


0.76 


1.5 0 


0 


0 


0.5 


0 


0 


0 


1 


0 


-0.089 -0.533 


0 


0 


0,22 


0.28 


0.6 1.5 


0 


0 


0.4 


0.5 


0 


0 


0 


1 


-0.044 -0.267 


0 


0 


0.4 


0 


0 0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1.333 0 


0 


0 


0.4 


0.4 


0 0 


0 


-1 


0 


0 


0 


0 


0 


0 


0.444 -1.333 


0 


0 


0.06 


-0.3 


-0.3 0 


0 


0' 


;p:,Z5 


0 


0 


0 


0 


0 


-0.089 -0.333 


1 


0 


0.03 


-0.04 


-0.4 -0.3 


0 


0 


0.25 


0.75 


0 


0 


0 


0 


0 -0.2 


0 


1 


Multiply Equation-7 by (-1/0.75) 


1.7 


0 


0 0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 0 


0 


0 


0.96 


1.7 


0 0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.311 -0.667 


0 


0 


0.44 


0.76 


1.5 0 


0 


0 


0.5 


0 


0 


0 


1 


0 


-0.089 -0.533 


0 


0 


0.22 


0.28 


0.6 1.5 


0 


0 


0.4 


0.5 


0 


0 


0 


1 


-0.044 -0.267 


0 


0 


0.4 


0 


0 0 


-1 


0 


0 


d 


0 


0 


0 


0 


-1.333 0 


0 


0 


0.4 


0.4 


0 0 


0 


-1 


0 


0 


0 


0 


0 


0 


0.444 -1 .333 


0 


0 


-0.08 


0.4 


0.4 0 


0 


0 


■■•■■-.1 


0 


0 


0 


0 


0 


0.1190.4444 


-1.333 


0 


0.03 


-0.04 


-0.4 -0.3 


0 


0 


0.25 


0.75 


0 


0 


0 


0 


0 -0.2 


0 


1 



Multiply Equation-5 by 0.5, add it to Equation-3 and replace Equation-3 

Page 15 



Multiply Equation-5 by 0.4. add it to Equation-4 and replace Equation-4 
Multiply Equation-5 by 0.25, add it to Equation-8 and replace Equation-8 



1.7 


0 


0 0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 


0 0 


0 


0.96 


1.7 


0 0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.311 


-0.667 0 


0 


0.4 


0.96 


1.7 0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.030 


-0.311 -0.667 


0 


0.188 


0.44 


0.76 1.5 


0 


0 


0 


0.5 


0 


0 


0 


1 


0.003 


-0.089 -0.533 


0 


0.4 


0 


0 0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1.333 


0 0 


0 


0.4 


0.4 


0 0 


0 


-1 


0 


0 


0 


0 


0 


0 


0.444 


-1 .333 0 


0 


-0.08 


0.4 


0.4 0 


0 


0 


-1 


0 


0 


0 


0 


0 


0.119 


0.444 -1.333 


0 


0.01 


0.06 


-0.3 -0.3 


0 


0 


0;-' 




0 


0 


0 


0 


0.030 


-0.089 -0.333 


1 



Multiply Equation-8 by (-1/0.75) 



1.7 


0 


0 0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 


0 


0 


0 


0.96 


1.7 


0 0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.311 


-0.667 


0 


0 


0.4 


0.96 


1.7 0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.030 


-0.311 


-0.667 


0 


0.188 


0.44 


0.76 1.5 


0 


0 


0 


0.5 


0 


0 


0 


1 


0.003 


-0.089 


-0.533 


0 


0.4 


0 


0 0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1.333 


0 


0 


0 


0.4 


0.4 


0 0 


0 


-1 


0 


0 


0 


0 


0 


0 


0.444 


-1.333 


0 


0 


-0.08 


0.4 


0.4 0 


0 


0 


-1 


0 


0 


0 


0 


0 


0.119 


0.444 


-1.333 


0 


-0.013 


-0.08 


0.4 0.4 


0 


0 


Of. 


;;j-i 


0 


0 


0 


0 


-0.040 


0.119 


0.444 


-1.333 


Multiply Equation-5 by 0.5, add it to Equation-4 and replace Equation-4 


1.7 


0 


0 0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 


0 


0 


0 


0.96 


1.7 


0 0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.311 


-0.667 


0 


0 


0.4 


0.96 


1.7 0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.030 


-0.311 


-0.667 


0 


0.181 


0.4 


0.96 1.7 


0 


0 


0 


0 


0 


0 


0 


1 


-0.017 


-0.030 


-0.31 1 


-0.667 


0.4 


0 


0 0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1.333 


0 


0 


0 


0.4 


0.4 


0 0 


0 


-1 


0 


0 


0 


0 


0 


0 


0.444 


-1.333 


0 


0 


-0.08 


0.4 


0.4 0 


0 


0 


-1 


0 


0 


0 


0 


0 


0.119 


0.444 


-1.333 


0 


-0.013 


-0.08 


0.4 0.4 


0 


0 


0 


-1 


0 


0 


0 


0 


-0.040 


0.119 


0.444 


-1.333 



Rearrange Equations 



1.7 


0 


0 


0 


0.667 0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


0.96 


1.7 


0 


0 


0.311 0.667 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0.4 


0.96 


1.7 


0 


0.030 0.311 


0.667 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0.181 


0.4 


0.96 


1.7 


0.017 0.030 


0.311 


0.667 


0 


0 


0 


1 


0 


0 


0 


0 


0.4 


0 


0 


0 


1.333 0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0.4 


0.4 


0 


0 


-0.444 1.333 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.08 


0.4 


0.4 


0 


-0.119 -0.444 


1.333 


0 


0 


0 


0 


0 


0 


0 


1 


0 



The new model coefficients with the PID dynamics removed are as follows: 
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Dependent Var-1 
Independent Var-1 Independent Var-2 

4,,= 0.96 4.2,2= 0.311 

4,3= 0.4 4.2,3= 0.030 

\,^,= 0.181 4,2,4= 0.017 

Dependent Var-2 
Independent Var-1 Independent Var-2 

4,u= 4,,,= 1.333 

4,1,2= 0-4 4,2,2= -0-444 

4,1,3= -008 4,2,3= -0.119 

4,1,4= -0.0133 4,2,4= 0.040 

[00043] Note that all the coefficient values changed. This new controller now has the dynamics of 
the second independent variable (a PID set point) removed. This controller can now be used to 
control the process and the development of this controller was done off line without having to do an 
additional time consuming expensive identification test on the process. 



Algorithm to Remove PID Dynamics, Open-loop Step Response fonm 

[00044] In the derivation and example, we discussed the algorithm to remove PID dynamics from an 
FIR model based on the impulse, or derivative, form of the equations. A similar algorithm can be 
derived for the Step coefficient form of the model, SO = AM , as well, as we will now illustrate with 
a 2 independent, 2 dependent variable example problem. For purposes of this example, we will 
assume that the second independent and second dependents are to be exchanged. The problem 
can be written in matrix notSition as: 



^1,1 








'lO' 






X 






X 


^2,1 ^2,2 








01 





[00045] We perform elimination operations (pivoting) to get: 



'Ai 0 " 




"a/V 




0-4.2 




'so,' 




X 






X 








A/, 
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[00046] Rearranging, we get: 







"a7," 




7o' 






X 








X 




_^2,1 ^2,2 _ 








_0/_ 




A/2 _ 



[00047] Which can be written as: 

[00048] Recall that the Impulse coefficients are defined as: 

l>ij,k = ^.,j,k - ^iMk-i) = ^/ j,i for ^ : 2 ^ ncoef 

[00049] Lil<ewise, we define the second difference coefficients as: 

^i.M=^M forA: = l 

= Kjjc - Kj,ik-i) = ^Kjj. for ^ : 2 -» ncoef 

[00050] Note that: 
Kj,m =Zky^) 



[00051] Note that the matrix is now a mixed bag of Step response coefficients {A), Impulse 
coefficients (6), and 2™^ difference coefficients (C). This is due to the fact that our new independent 
variable is in the "accumulative" form instead of the "delta" form and the new dependent variable is 
in the "delta" forni instead of the "accumulative" form. To convert this system of equation to the 
Step form and thereby recover the Step coefficients, two steps must be performed: 



a 
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[00052] Step-1 : Convert new independent variable from "accumulative" to "delta" form, 

[00053] Step-2: Convert new dependent variable from "delta" to "accumulative" form, nJ^^dl^ 

[00054] Step-1 : Convert new independent variable from "accumulative" to "delta" form. 

[00055] This step requires only a rearrangement of terms in the equations. Note that<502 appears 
two sections of the matrix: 




^2,0) 

<52,o)+ 4,2.1 (<^2.2 
<^2.o)+ 4,2,2 fe,2 



^^2,0 )+ 4,2.1 (<52,3 -^2,0) 




02,0) 

^2,0 ^2,2,1 (^2,3 ~ ^2,0 ) 



Since 




forA: = l 



= Aa^.yjt for kil-^ncoef 



and 



fork = l 



= Ab^jf^ fork:2^ncoef 



we can write the above as: 
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^1,2,1 (<^2.l""<^2.o) 
(«1,2.2 - «1.2.1 X<^2,1 - <^2,0 ) + «U2,1 (^2,2 " <^2,0 ) 

(^1,2,3 - «1,2,2 fej - <^2.0 K (^1,2,2 " '^1,2.1 X^2,2 " <^2.0 )+ ^1,2,1 (^2,3 " <^2,0 ) 

4,2,1(^2.1 -<^2.o) 
fe,2.2- 4,2,1 fo,l-C>2,o) +4,2,1^2-^^2,0) 

(4.2,3 -4.2,2^,1 -<^2.oK(4.2.2 -4.2.l)(C>2.2 " <^2,0 K 4.2.1 (^2.3 -<^2,o) 



«I,2,l(<^2,|-<^2,o) 

^1,2,2<^2,1 -^1,2,2<^2,0 -^1,2,1<^2,1 +^l,2,l'^2,C +^1,2,1<^2,2 -'2l,2,l<^2,0 

^1,2,3^2,1 -Klpifi -^1.2,2^2,1 +«I.2.2<52.0 +^1,2,2<52,2 -«1,2,2^2,0 -«I,2,1<52,2 +«!,2,1<52,0 +^1,2,1<^2,3 -«1,2.1<^2.0 
4,2,1 ('^2,1 -^2,0) 

4,2,2*^2,1 -\2.2<52,0 -^2,2,I<52,1 +^2,2,1<52,0 +^2,2,1<^2,2 -^2,2,1<^2,0 

4,2.3C>2.1 -kxi^2,0 -4,2.2(^2,1 +4.2,2C>2,0 +4,2,2C'2.2 -4.2,2^*2,0 'KlPia ^KlPlfi +4,2,1^^2,3 -^2.2,|C'2.0 
«1.2.l(C>2,l-<52,o) 

«1,2,2(02., - <52,o)+ «l,2,l(02,2 " ^^^.l + O2.O " <52,o) 

«1,2,3(C'2,1 - ^2,0)+ «I.2,2{C>2,2 " C>2,1 + ^2,0 " <^2,o)+ «1.2,l{<52,3 " ^2.2 + <?2,0 " C>2,o) 

^2,2,1(^2,1 ~ ^2,0) 

4,2,2(02.1 -(52,o)+4,2,l(C>2,2 -<52,I +^2,0 -<52,o) 

4,2,3(02,1 - 02,0)+ 4,2.2(02,2 - O2,, + 02,0 - 02,0)+ 4,2,1(02,3 - O2.2 + 0^,0 - 



4,2^^2 
C2 2^ O2 



«l,2,l(02,I -02,0) 

«i,2.2 (O2,, - 02,0 ) + a,,2,i (02,2 - O2,, ) 

«1,2,3 (02,1 - O2.0 ) + ^1,2,2 (O2.2 - O2,, ) + a, 2,1 (0^,3 - ^ ) 

4,2,1(02,1 —Ojo) 

4,2.2 (O2,, - 02_o ) + 4.2,1 (02,2 - O2,, ) 

4.2.3(02.1 - O2.0) + 4,2.2(02.2 - O2,,) + 4.2.1(02.3 - 02,2) 
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«U,2^^2.1+ ^1.2.1^^^2,2 

«1,2,3^^^2,1 + <3l.2,2^^^2,2 + ^1,2,1^^^2.3 

^2.2,l'^^^2,l 

4,2,2AC>2,1 +4.2.1^^^2.2 

4.2,3^C>2,1 + kxi^^ia + 4,2,1^^^2.3 



A,2A0, 
4,2AO, 



[00056] Since 



, we can re-write the system of equations as: 



[00057] This completes Step-1 . 

[00058] Step-2: Convert new dependent variable from "delta" to "accumulative" fomi. 



[00059] The equations for the new second dependent variable are written below. It is necessary to 
accumulate these equations to convert from the "delta" to the "accumulative" form. 



4,U^U+4,lA+4lA3 +4.2,3^4+4,2,2^2+4,2,1^3 
4iAi^-4,1,3A^1,2+4A+4,U^M 



=4l-4o=^2,l 
=^2,2~"^,1^^2,2 
=72^3-/2^=^/23 



Since by definition, 6,^^ j = a^^..^ , and i -/^^o = ^7,1 = ^ ^7,1 ' ^^^^ ^'^^^ equation becomes; 

^2,1,1^1,1+^2,2,1^^2,1 =<^^2,l 
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[000601 To obtain the second Step coefficient equation, add together the first two Impulse coefficient 
equations: 

(4,1,1 +4,1,2)^1,1 +4,1,1^1,2 +(4,2,1 +4,2,2)^02,1 +4.2,1^^2.2 =^2.2 -^,1 +^2,1 =^2.2 "^2,0 

or, 

^2,1,2^1,1 +^2,M^l,2 +'^2,2,2^<^2,1 +^2.2,1^^2.2 =^^2,2 

[00061] To obtain the third Step coefficient equation, add together the first three Impulse coefficient 
equations: 

(4.1,1 +4,1.2 +4,1,3)^1,1 +(4.1,1 +4,1,2)^1,2 +4.1,lA/i,3 

+ (4,2,1 +4.2.2 +4,2,3 )^<^2.l +(4,2,1 +4,2,2)^02.2 +4,2.1^02,3 

— -^2.3 ~^2,2 +-^2,2 ~^2,1 +"^2.1 ~^2,0 ~ ^2,3 ""-^2,0 

or, 

^2,1.3^1,1 +^2,1.2^1.2 +^2,1,1^1,3 +«2,2.3^02,l +^2,2,2^02,2 +«2,2,1^023 = ^ ^z.S 

[00062] To obtain the fourth Step coefficient equation, add together the first four Impulse coefficient 
equations: 

(4,1.1 +4.1.2 +4,1,3 +4.1,4)^1,1 +(4.1.1 +4.1.2 +4.1.3)^1.2 +(4.1,1 +4,1,2)^1.3 +4.1,lA/i,4 

+ (4.2.1 +4.2.2 +4.2.3 +4,2.4)^02,, +(4,2.1 +4.2.2 +4,2.3)^02,2 +(4.2,1 +4,2,2)^02,3 +4,2.1^02,4 

= ^2,4 ~^2.3 +-^2,3 "^2.2 +^2,2 "^2.1 +^2,1 ~^2.0 ="^2,4 ~-^2.0 

or, 

fl2,j,4A/i, +^2,1.3^1,2 +^2,1,2^1.3 +^2,1,1^1,4 +^2,2,4^02,i +^2.2.3^02,2 +^2,2,2^02,3 +^2,2,1^02,4 ==SI^^ 

[00063] So the system of equations for the new 2""^ dependent variable now becomes: 

^2,UlM,l +^2,2,1^0, 1 = ^2,1 ~ 4,0 = ^ h^ 

^2,1,2^1,1 + ^2,l.lM,2 + ^2^.2^^2,1 + ^1,2^^2,2 " 4.2 " 4.0= ^ 4,2 

^2,1,3^,1 + ^2,1,2^/1,2 + ^2,l,lM.3 + ^2,2,3^02,1 + €1,^2^0^,1 + ^2,2,1^03,3 = 4,3 " 4,o = S I^ ^ 

£2,1,4 A/1,1 + ^22,1.3^1,2 + ^2.1,2^1,3 + ^2,1,1^1,4 + ^2,2,4^02,, + a^x^^^^.^^ + ^2,2.2^^2,3 + CL^xi^O^A = 4,4 ^ 4.0 = ^ 4.4 
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[00064] And the overall system of equations becomes: 



• • • • 


4^^^! -fV^ 

• • • • 

• • • • 


— 

=^ 


^ AT 

• • • • 


+a^;y/!qi -^1^ 


-^2,1 



i3 which can be re-written as: 



[00065] To further illustrate the application of this invention another numerical example is given to 
demonstrate the use of the algorithm just derived for the open-loop step response model. This 
algorithm is applied to equations of the form SO = AM . Given a model with two (2) independent 
variables, two (2) dependent variables and four (4) model coefficients, where the second 
independent variable is the set point of a PID controller and the second dependent variable is the 
valve position of the PID controller, we wish to re-constitute the model with the PID valve position 
as an independent variable instead of the PID set point. This requires that the dynamics of the PID 
controller be removed from all system responses according to the algorithm previously discussed. 
The underlying model in this example is the same as that used in Appendix-2. 



Dependent Var-1 

Independent Var-1 Independent Var-2 



^1,1,1 ~ 


1.5 


^1.2,1 


0.5 


^\,\,2 ~ 


2.1 


'^1,2,2 ~ 


0.9 


%3 = 


2.3 


'^1,2,3 ~ 


1.1 


^1,1,4 = 


2.4 


'^1,2,4 ~ 


1.2 
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Dependent Var-2 

Independent Var-1 Independent Var-2 



^2,1,1 ~ 


-0.3 


^2.2,1= 0.75 


^2,1,2 ~ 


-0.7 


flj 2 2 ~ ^ -0 


^2,1,3 ~ 


-0.8 


«2,2,3~ '^•''^ 


^2,1,4 ~ 


-0.85 


«2,2,4=1-2 



[000661 The problem is specified in tlie matrix below. 



' Indicates Pivot Element 



1.5 


0 


0 


0 


0.5 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


2.1 


1.5 


0 


0 


0.9 


0.5 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


2.3 


2.1 


1.5 


0 


1.1 


0.9 


0.5 


0 


0 


0 


1 


0 


0 


0 


0 


0 


2.4 


2.3 


2.1 


1.5 


1.2 


1.1 


0.9 


0.5 


0 


0 


0 


1 


0 


0 


0 


0 


-0.3 


0 


0 


0 


;0.75 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.7 


-0.3 


0 


0 


1 


0.75 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.8 


-0.7 


-0.3 


0 


1.15 


1 


0.75 


0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.85 


-0.8 


-0.7 


-0.3 


1.2 


1.15 


1 


0.75 


0 


0 


0 


0 


0 


0 


0 


1 


Multiply Equation-5 by (-1/0.75) 


1.5 


0 


0 


0 


0.5 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


2.1 


1.5 


0 


0 


0.9 


0.5 


0 


d 


0 


1 


0 


0 


0 


0 


0 


0 


2.3 


2.1 


1.5 


0 


1.1 


0.9 


0.5 


0 


0 


0 


1 


0 


0 


0 


0 


0 


2.4 


2.3 


2.1 


1.5 


1.2 


1.1 


0.9 


0.5 


0 


0 


0 


1 


0 


0 


0 


0 


0.4 


0 


0 


0 




0 


0 


0 


0 


0 


0 


0 


-1.333 


0 


0 


0 


-0.7 


-0.3 


0 


0 


1 


0.75 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.8 


-0.7 


-0.3 


0 


1.15 


1 


0.75 


0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.85 


-0.8 


-0.7 


-0.3 


1.2 


1.15 


1 


0.75 


0 


0 


0 


0 


0 


0 


0 


1 



Multiply Equation-5 by 0.5, add it to Equation-1 and replace Equation-1 

Multiply Equation-5 by 0.9, add it to Equation-2 and replace Equation-2 

Multiply Equation-5 by 1.1, add it to Equation-3 and replace Equation-3 

Multiply Equation-5 by 1 .2, add it to Equation-4 and replace Equation-4 

Multiply Equation-5 by 1 .0, add it to Equation-6 and replace Equation-6 

Multiply Equation-5 by 1.15, add it to Equation-7 and replace Equation-7 

Multiply Equation-5 by 1.2, add it to Equation-8 and replace Equation-8 



1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 


0 


0 


0 


2.46 


1.5 


0 


0 


0 


0.5 


0 


0 


0 


1 


0 


0 


-1.200 


0 


0 


0 


2.74 


2.1 


1.5 


0 


0 


0.9 


0.5 


0 


0 


0 


1 


0 


-1.467 


0 


0 


0 
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2 1 


1 5 


0 


1.1 


0.9 


0.5 


0 


0 


0 


1 


-1 .600 


0 


0 


0 




u 


u 




1 


n 


n 




n 

V 


0 


0 


0 


-1 333 

i .SJ\J\J 


0 


0 


0 


-0.3 


-0.3 


0 


0 


0 


0.75 


0 


0 


0 


0 


0 


0 


-1.333 


1 


0 


d 


-0.34 


-0.7 


-0.3 


0 


0 


1 


0.75 


0 


0 


0 


0 


0 


-1.5 


0 


1 


0 


-0.37 


-0.8 


-0.7 


-0.3 


0 


1.15 


1 


0.75 


0 


0 


0 


0 


-1.600 


0 


0 


1 


Multiply Equation-6 by (-1/0.75) 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 


0 


0 


0 


2.46 


1.5 


0 


0 


0 


0.5 


0 


0 


0 


1 


0 


0 


-1.200 


0 


0 


0 


2.74 


2.1 


1.5 


0 


0 


0.9 


0.5 


0 


0 


0 


1 


0 


-1.467 


0 


0 


0 


2.88 


2.3 


2.1 


1.5 


0 


1.1 


0.9 


0.5 


0 


0 


0 


1 


-1 .600 


0 


0 


0 


0.4 


0 


0 


0 


-1 


0 


0 






0 


0 


0 


-1.333 


0 


0 


0 


0.4 


0.4 


0 


0 


0 MM-.-"'' 


0 


: 




0 


0 


0 


1.778 


-1.333 


0 


0 


-0.34 


-0.7 


-0.3 


0 


0 


1 


0.75 




J 


0 


0 


0 


-1.533 


0 


1 


0 


-0.37 


-0.8 


-0.7 


-0.3 


0 


1.15 


1 


0.75 




0 


0 


0 


-1.600 


0 


0 


1 



Multiply Equation-5 by 0.5, add it to Equation-2 and replace Equation-2 
Multiply Equation-5 by 0.9, add It to Equation-3 and replace Equation-3 
Multiply Equation-5 by 1.1, add it to Equation-4 and replace Equation-4 
Multiply Equation-5 by 1.0, add it to Equation-7 and replace Equation-7 
Multiply Equation-5 by 1.15, add it to Equation-8 and replace Equation-8 



1.7 0 0 0 
2.66 1.7 0 0 

3.1 2.46 1.5 0 
3.32 2.74 2.1 1.5 


0 0 0 0 
0 0 0 0 
0 0 0.5 0 
0 0 0.9 0.5 


10 0 0 
0 10 0 
0 0 10 
0 0 0 1 


-0.667 0 0 0 
-0.311 -0.667 0 0 
0.133-1.200 0 0 
0.356 -1 .467 0 0 


0.4 0 0 0 
0.4 0.4 0 0 

0.06 -0.3 -0.3 0 
0.09 -0.34 -0.7 -0.3 


-10 0 0 
0-100 

0 0 "OjfS 0 
0 0 1 0.75 


0 0 0 0 

0 0 0 0 
0 0 0 0 
0 0 0 0 


-1.333 0 0 0 

1.778-1.333 0 0 
0.244-1.333 1 0 
0.444-1.533 0 1 


Multiply Equation-7 by (-1/0.75) 


1.7 0 0 0 
2.66 1.7 0 0 

3.1 2.46 1.5 0 
3.32 2.74 2.1 1.5 


0 0 0 d 
0 0 0 0 

0 0 0.5 0 
0 0 0.9 0.5 


10 0 0 
0 10 0 

0 0 10 
0 0 0 1 


-0.667 0 0 0 
-0.311 -0.667 0 0 
0.133-1.200 0 0 
0.356 -1 .467 0 0 


0.4 0 0 0 
0.4 0.4 0 0 
-0.08 0.4 0.4 0 
0.09 -0.34 -0.7 -0.3 


-10 0 0 
0-100 

0 0 -1 o 

0 0 1 0.75 


0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 


-1.333 0 0 0 
1.778-1.333 0 0 

-0.326 1.778 -1.333 0 
0.444 -1.533 0 1 



Multiply Equation-5 by 0.5, add it to Equation-3 and replace Equation-3 
Multiply Equation-5 by 0.9, add it to Equation-4 and replace Equation-4 
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Multiply Equation-5 by 1 .0, add it to Equation-8 and replace Equation-8 



1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 0 0 


0 


2.66 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.311 -0.667 0 


0 


3.06 


2.66 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.030 -0.311 -0.667 


0 


3.248 


3.1 


2.46 


1.5 


0 


0 


0 


0.5 


0 


0 


0 


1 


0.062 0.133-1.200 


0 


0.4 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1 .333 0 0 


0 


0.4 


0.4 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


0 


1.778-1.333 0 


0 


-0.08 


0.4 


0.4 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


-0.326 1.778-1.333 


0 


0.01 


0.06 


-0.3 


-0.3 


0 


0 


0-: 


'd>5 


0 


0 


0 


0 


0.119 0.244-1.333 


1 


Multiply Equatlon-8 by (-1/0.75) 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 0 0 


0 


2.66 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.311 -0.667 0 


0 


3.06 


2.66 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.030 -0.311 -0.667 


0 


3.248 


3.1 


2.46 


1.5 


0 


0 


0 


0.5 


0 


0 


0 


1 


0.062 0.133-1.200 


0 


0.4 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1.333 0 0 


0 


0.4 


0.4 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


0 


1 .778 -1 .333 0 


0 


-0.08 


0.4 


0.4 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


-0.326 1.778-1.333 


0 


-0.013 


-0.08 


0.4 


0.4 


0 


0 


0 


■■■■ ,-1 


0 


0 


0 


0 


-0.158-0.326 1.778 


-1.333 


Multiply Equation-5 by 0.5, add it to Equation-4 and replace Equation-4 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


-0.667 0 0 


0 


2.66 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.311 -0.667 0 


0 


3.06 


2.66 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.030 -0.311 -0.667 


0 


3.241 


3.06 


2.66 


1.7 


0 


0 


0 


0 


0 


0 


0 


1 


-0.017 -0.030 -0.311 


-0.667 


0.4 


0 


0 


0 


-1 


0 


0 


0 


0 


0 


0 


0 


-1 .333 0 0 


0 


0.4 


0.4 


0 


G 


0 


-1 


0 


0 


0 


0 


0 


0 


1.778-1.333 0 


0 


-0.08 


0.4 


0.4 


C 


0 


0 


-1 


c 




0 


0 


c 


-0.326 1.778-1.333 


0 


-0.013 


-0.08 


0.4 


0.4 


0 


0 


0 




0 


0 


0 


c 


-0.158-0.326 1.778 


-1.333 



Rearrange Equations 



1.7 


0 


0 


0 


0.667 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


0 


2.66 


1.7 


0 


0 


0.311 


0.667 


0 


0 


0 


1 


0 


0 


0 


0 


0 


0 


3.06 


2.66 


1.7 


0 


0.030 


0.311 


0.667 


0 


0 


0 


1 


0 


0 


0 


0 


0 


3.241 


3.06 


2.66 


1.7 


0.017 


0.030 


0.311 


0.667 


0 


0 


0 


1 


0 


0 


0 


0 


0.4 


0 


0 


0 


1.333 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 


0.4 


0.4 


0 


0 


-1.778 


1.333 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


-0.08 


0.4 


0.4 


0 


0.326 


-1.778 


1.333 


0 


0 


0 


0 


0 


0 


0 


1 


0 


-0.013 


-0.08 


0.4 


0.4 


0.158 


0.326 


-1.778 


1.333 


0 


0 


0 


0 


0 


0 


0 


1 
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Accumulate coefficients for new 2""^ independent variable 



1.700 0.000 0.000 0.000 
2.660 1.700 0.000 0.000 

3.241 3.060 2.660 1.700 


0.667 0.000 0.000 0.000 
0.978 0.667 0.000 0.000 

1 007 0 Q78 0 667 0 000 
1.024 1.007 0.978 0.667 


0.400 0.000 0.000 0.000 
0.400 0.400 0.000 0.000 
-0.080 0.400 0.400 0.000 
-0.013 -0.080 0.400 0.400 


1.333 0.000 0.000 0.000 
-0.444 1.333 0.000 0.000 
-0.119 -0.444 1.333 0.000 

0.040 -0.119-0.444 1.333 



Accumulate coefficients for new 2 independent variables 



1.700 0.000 0.000 0.000 
2.660 1.700 0.000 0.000 

3.060 2.660 1.700 0.000 
3.241 3.060 2.660 1.700 


0.667 0.000 0.000 0.000 
0.978 0.667 0.000 0.000 
1.007 0.978 0.667 0.000 
1.024 1.007 0.978 0.667 


0.400 0.000 0.000 0.000 
0.800 0.400 0.000 0.000 
0.720 0.800 0.400 0.000 
0.707 0.720 0.800 0.400 


1.333 0.000 0.000 0.000 
0.889 1.333 0.000 0.000 

0.770 0.889 1.333 0.000 
0.810 0.770 0.889 1.333 



[00067] The new model coefficients with the PID dynamics removed are as follows: 

Dependent Var-1 



Independent Var-1 



= 1.700 



1 2 = 2.660 
^113= 3.060 
^11,4= 3.241 



Dependent Var-2 

Independent Var-1 

a^,, =0.400 

^2,1,2=0-800 
^2,1,3=0.720 
«2,i,4 =0-707 

Note that all the coefficient values changed. 



Independent Var-2 
a^2^= 0.667 

= 0.978 



a, 



1,2,2 



ai,2,4= 1-024 



Independent Var-2 
=0.889 



a. 



2,2,2 



«2.2,3 =0-770 
a2 2,4 =0.810 



Check that corresponding Impulse coefficients match those identified with the FIR example. 



Page 27 



Dependent Var-1 



Independent Var-1 



Independent Var-2 



^1,1,2= 0.960 
Z?ii3= 0.400 



\^,= 0.667 

^1,2,2= 0.311 

^1,2,3= 0-030 
^2,4= 0.017 



Dependent Var-2 



Independent Var-1 



Independent Var-2 



^2,1,1= 0.400 
^2,u = 0.400 
624,3= -0.080 
^12,1,4= -0013 



^12,24- 1-333 
^^2,2,2= -0.444 

^2,3= -0.119 

^2,2,4= 0.040 



Column Simulation Example 



[00068] Yet another embodiment of the use of the algorithm is demonstrated in the following 
example. This example will illustrate the following: 

[00069] The use of a valve-based Finite Impulse Response (FIR) model as a process simulator. 
[00070] Plant step-test and Identification of an FIR model based on a specific regulatory control 
configuration. 

[00071] Use of the proposed algorithm to remove the PID controller dynamics and recover the 
underlying valve-based model. 

[00072] In this example, an FIR model based on valve positions is used as the process model to 
simulate the behavior of a complex fractionator. The regulatory control for the fractionator consists 
of three PI (proportional/integral) feedback controllers. A plant step test is performed on the 
simulation using the regulatory controller set points. An FIR model is then obtained for the 
fractionator based on the set points of the PI controllers. This model based on the regulatory 
control scheme is then input to the algorithm to remove the PI controller dynamics and recover the 
original FIR process model. 

[000731 It should be noted that the term Finite Impulse Response (FIR) model is used to refer to the 
open-loop step response form of the models, since the step form could be directly calculated from 
the impulse coefficients. 
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[00074] Description of Complex Fractionator Schematic 

[00075] The schematic for the Complex Fractionator is shown in Figure 2. The feed flow rate 5 is 
controlled by the upstream unit and is pre-heated in a furnace 6. The fractionator 7 has a top, 
middle and bottom product. The fractionator overhead temperature is controlled with a PI controller 
8 moving the top reflux. The middle product draw temperature is controlled with a PI controller 9 
moving the middle product draw rate. A third PI controller 10 moves the bottom product rate to 
control the fractionator bottoms level. The bottom composition (light component) is measured with 
an analyzer 11. 

[00076] Description of Finite Impulse Response (FIR) Model 

[00077] The process model used in this example is an open-loop, step response model based on the 
valve positions, summarized as follows: 

Model Independent Variables 
TIC-2001 .OP - Top Reflux Flow Valve 
TIC-2002.OP - Middle Product Flow Valve 
LIC-2007.OP - Bottoms Product Flow Valve 
FIC-2004.SP - Middle Reflux Flow Rate 
FI-2005.PV - Fractionator Feed Rate 

Model Dependent Variables 
TIC-2001 .PV - Fractionator Overhead Temperature 
TIC-2002.PV - Middle Product Draw Temperature 
LIC-2007.PV - Fractionator Bottoms Level 

AI-2022.PV - Fractionator Bottoms Composition (Light Component) 

[00078] The open-loop step response model can be viewed in an idealized sense as being 
generated as follows. With the system at steady state, the first Independent variables is increased 
by one engineering unit at time=0 while holding all other independent variables constant. The 
values for all dependent variables are then measured at equally spaced time intervals until the 
system reaches steady state again. The model response curves for each dependent variable with 
respect to the first independent variable are then calculated by subtracting the value of the 
dependent variable at time=0 from each of the measured values at each future time interval for that 
dependent variables. Essentially, a step response curve represents the effect on the dependent 
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variable of a change in the independent variable. This process is then repeated successively for all 
the independent variables to generate the full model. The steady state time for the model is defined 
by the steady state time of the slowest response curve in the system. 

[000791 Clearly in the real world, the model cannot be generated in this fashion since often the 
process is not at steady state. Further, it is impossible to prevent measured and unmeasured 
disturbances from affecting the system during an independent variable step. Generation of the 
model requires that multiple steps be made in each independent variable (plant step test). The 
data thus collected is then analyzed with a software package such as AspenTech's DMCplus Model 
program to calculate the open-loop step response model. 

[00080] Once such a model has been identified, it can be used to predict future system response 
based on past changes in the independent variables. That is to say, if we know how all 
independent variables have changed for one steady-state time into the past, we can use the model 
to predict how the dependent variables will change for one steady-state time into the future, 
assuming no further independent variable changes. This illustrates the use of the model for 
Prediction. (This is the basis for using an FIR model as a process simulator). 

[00081] Given the predicted future system response based on no further independent variable 
changes and given the constraints on all independent and dependent variables, the model can be 
used to plan a strategy of independent variable moves to keep all independent and dependent 
variables within constraints. This illustrates the use of the model for Control. 

Using a Finite Impulse Response (FIR) Model as a Process Simulator 

[00082] The model for this example has a steady state time of ninety (90) minutes. A three (3) 
minute time interval is used. The resulting response curves are each defined by a vector of thirty 
(30) numbers representing the accumulative change in that dependent variable across time with 
respect to a step change in the independent variable at time=0 while holding all other independent 
variables constant. 

[00083] The model coefficients are shown in Table 1 and the model plots are shown in Figure 3. 
This model, based on valve positions, is used to predict future system behavior in the model 
dependent variables based on past and present changes in the model independent variables. 
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Table 1: Fractionator Simulation Valve-based Model Coefficients 



Step Response Coefficients for Dependent Variable-1 : TIC-2001 .PV DEG F 





TIC-2001. OP 


TIC-2002.OP 


L1C-2007.OP 


FIC-2004.SP 


FI-2005.PV 




+1 % Move 


+1 % Move 


+1 % Move 


+1 MBBL/D Move 


+1 MBBUD Move 


Minutes 


at Time=0 


at Time=0 


at Time=0 


at Time=0 


at Time=0 


0 


0.000 


0.000 


0.0 


0.00 


00 


3 


-0.101 


-0.048 


0.0 


-2.05 


2.9 


6 


-0.175 


-0.076 


0.0 


-3.58 


6.1 


9 


-0.206 


-0.088 


0.0 


-4.43 


7.5 


12 


-0.227 


-0.068 


0.0 


-5.03 


7.8 


15 


-0.245 


-0.040 


0.0 


-5.58 


8.2 


18 


-0.262 


-0.015 


0.0 


-6.16 


8.5 


21 


-0.277 


0.010 


0.0 


-6.65 


8.6 


24 


-0.292 


0.033 


00 


-7 04 


89 


27 


-0.306 


0.054 


00 


-7.37 


90 


30 


-0.323 


0.069 


00 


-7.67 


9.3 


33 


-0.340 


0.084 


0.0 


-7.95 


9.5 


36 


-0.356 


0.096 


0.0 


-8.18 


96 


39 


-0.372 


0105 


0.0 


-8.37 


9.8 


42 


-0.386 


0.113 


0.0 


-8.52 


9.8 


45 


-0.399 


0.121 


0-0 


-6.65 


9-8 


48 


-0.410 


0.128 


0.0 


-8.75 


9.9 


51 


-0.420 


0.135 


0.0 


-8 84 


10.0 


54 


-0.428 


0.140 


00 


^92 


10.1 


57 


-0.435 


0.145 


0.0 


-^.98 


10.3 


60 


-0.440 


0.149 


0.0 


-9.04 


10.4 


63 


-0.445 


0.153 


0.0 


-9.09 


10.5 


66 


-0 450 


0,156 


0.0 


-9.13 


10.5 


69 


-0.453 


0159 


0.0 


-9.17 


10.5 


72 


-0.457 


0 161 


0.0 


-9.21 


10.5 


75 


-0.460 


0.163 


0.0 


-9.24 


10.4 


78 


-0.462 


0.165 


0.0 


-9.26 


10.4 


81 


-0.464 


0.166 


0.0 


-9.28 


10.4 


84 


-0.465 


0.167 


0.0 


-9.29 


10.4 


87 


-0.466 


0.167 


00 


-9.29 


10.4 


90 


-0.466 


0.167 


0.0 


-9.29 


10.5 


Step Response Coefficients for Dependent Variab 


e-2: TIC-2002.PV DEG F 




TIC-2001. OP 


TIC-2002.OP 


UC-2007.OP 


FIC-2004.SP 


FI-2005.PV 




+1 % Move 


+1 % Move 


+1 % Move 
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Step Response Coefficients for Dependent Variable-3: LIC-20Q1 .PV % 
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e-4: AI-2022PV MOLE % 
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[0084] As mentioned above, there are three PI (Proportional/Integral) controllers in the system. 
These PI controllers are configured as follows: 



Table 2: Fractionator PID Controllers 



PID Loop Name 


Set Point 


Process 
Variable 


Output 


Kp 


Ki 


Top Temperature 


TIC-2001.SP 


TIC-2001.PV 


TIC-2001.OP 


-2.0 


3.0 


Middle Product Draw 
Temperature 


TIC-2002.SP 


TIC-2002.PV 


TIC-2002.OP 


3.0 


8.0 


Bottoms Level 


LIC-2001.SP 


LIC-2001.PV 


LIC-2007.OP 


-1.0 


4.0 



[0085] A plant test was performed (data plots in FIG. 4) with these PI controllers regulating the 
process. The independent and dependent variables for the system were as follows: 



Model Independent Variables 

TIC-2001 .SP - Top Reflux Flow Valve SP 
TIC-2002.SP - Middle Product Flow Valve SP 
LIC-2007.SP - Bottoms Product Flow Valve SP 
FIC-2004.SP - Middle Reflux Flow Rate 
FI-2005.PV - Fractionator Feed Rate 



Model Dependent Variables 

Tl C-200 1 . PV - Fractionator Overhead Temperature 

TIC-2002.PV - Middle Product Draw Temperature 

LIC-2007.PV - Fractionator Bottoms Level 

TIC-2001 .OP - Top Reflux Flow Valve 

TIC-2002.OP - Middle Product Flow Valve 

LIC-2007.OP - Bottoms Product Flow Valve 

AI-2022.PV - Fractionator Bottoms Composition (Light Component) 

[0086] This illustrates the use of a valve-based FIR model as a process simulator. As described 
above, the PID control calculations were performed external to the process simulation. 
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[0087] The resulting data were analyzed and a model based on this PID configuration was 
identified, as shown in FIG. 5. 

[0088] The new algorithm to remove PID dynamics was applied to the model shown in FIG. 5, and 
this model with the PID dynamics removed is compared to the original simulation model. As can be 
seen in FIG. 6, the algorithm successfully recovers the original valve based model. Note that the 
steady state time of the recovered model is longer than the steady state time of the original model. 
This is a result of a longer steady state time for the model with the PID controllers. The original 
valve-based simulation model had a steady state time of 90 minutes. When the PID controllers 
were configured and the plant step-test performed, it took 180 minutes for the process to reach 
steady state, due to having to wait for the PID feedback control to settle out. The steady state time 
of the recovered valve-based model has the same steady state time as the model containing the 
PID dynamics from which it was generated. It can be seen, however, that the recovered model has 
reached steady state in 90 minutes, and if it were truncated at that point, would exactly match the 
original valve-based model. 
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Applicability of the Invention 

10089] In the past, when the PID controllers were re-tuned or when the regulatory control scheme 
was reconfigured, a new plant was performed and a new model constructed. The invention 
described in this document removes the PID controller dynamics without having to perform another 
plant test. 

[0090] This ability to remove PID dynamics allows creation of an off-line process simulator based 
only on valve positions instead of PID set points. The plant test can be performed with any stable 
regulatory configuration and PID tuning and a corresponding model can be obtained. The algorithm 
to remove the PID dynamics Is then applied to the resulting model to remove the dynamics of all 
PID controllers and convert the model inputs from set points to valves. The regulatory control 
scheme can then be emulated external to the process model via a DCS console or console 
emulator. This allows the operator to put PID controllers in Manual-mode, break cascades, retune 
PID controller, or even re-configure the regulatory control scheme. 

[0091] With regard to model-based control applications, there are times when it is necessary to 
modify the PID tuning of a PID controller in the system. With the ability to remove the PID 
dynamics, a model can be generated which is based on the valve of this PID controller. The off-line 
simulation calculation can then be performed to generate a new process model that contains the 
new PID tuning, and this updated model can be incorporated into the model based controller, thus 
preventing a plant step test. This technique can also be applied if the regulatory control scheme is 
to be reconfigured. Assume that we have a temperature controller set point as an input to our 
model. If that valve is stuck and cannot be repaired without shutting the unit down, the algorithm 
could be applied to remove the dynamics of the temperature controller and the control application 
could continue to be used without the temperature controller. 

[0092] Another advantage of this invention is that a process can be tested in one regulatory 
configuration and a model-based controller can be commissioned with a different configuration. An 
example is a Fluidized Bed Catalytic Cracking Unit (FCCU) where the system pressure is controlled 
with a PID controller moving the speed of the Wet Gas Compressor. Often the most economical 
place to run the unit is with the compressor at maximum speed, but in this case, the pressure is not 
directly controlled. Testing the unit with the pressure off control is difficult. The solution is to test 



Page 35 



the plant with the PID controller moving the compressor speed, keeping the speed on control. 
When the model is obtained, the pressure controller PID dynamics are removed and the model 
based control application will them move the compressor speed directly. In this example, the model 
based control application controls the system pressure as an output by manipulating other inputs 
when the compressor speed is at maximum. 

[0093] Often when testing a unit, valves of certain PID controllers are driven off control during the 
plant test. At the present time, this data cannot be used in constructing the process model. With 
the new algorithm, it is possible to use all the data, even when a PID controller is off control. This is 
done by first identifying the model as before using data only where the PID controller is on control. 
This model is then modified to remove the PID dynamics and the new data is "filtered into" the 
model. 

[0094] Thus, this new invention will allow construction of high fidelity, useable off-line process 
simulators and will enhance the ability to implement and maintain model-based control applications. 

[0095] While a preferred form of the invention has been disclosed and described in the drawings, 
since variations in the preferred form will be evident to those skilled in the art, the invention should 
not be construed as limited to the specific forms shown and described, but instead is as set forth in 
the following claims when read in the light of the foregoing disclosure. 
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